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ABSTRACT 

Disk accretion onto a slowly rotating neutron star with a weak magnetic field 
H < 3 X 10^ gauss is considered in a wide range of luminosities 1/100 < L/L^.^^ < 1; 
where L(,dd is the Eddington luminosity. We construct a theory for the deceleration 
of rotation and the spread of matter over the stellar surface in the shallow-water 
approximation. The rotation slows down due to friction against the dense underlying 
layers. The deceleration of Keplerian rotation and the energy release take place on 
the stellar surface in a latitudinal belt whose meridional width rises with increasing 
L. The combined effect of centrifugal force and radiation pressure gives rise to two 
latitudinal rings of enhanced brightness which are symmetric around the equator in 
the upper and lower hemispheres. They lie near the edges of differentially rotating 
and radiating upper and lower belts. The bright rings shift from the equatorial zone 
to higher latitudes when the luminosity L rises. The ring zones are characterized by a 
minimum surface density and, accordingly, by a maximum meridional spread velocity. 
At a low accretion rate and luminosity, the released energy is removed through the 
comptonization of low- frequency photons. 



1. OVERALL PICTURE 

In the Newtonian approximation, half of the gravitational energy is released in an extended 
disk during accretion onto a neutron star with a weak magnetic field. If the star rotates slowly, 
the other half of the energy must be released in the immediate vicinity of the star as the matter 
transfers from the disk onto the surface. This energy release is generally thought to occur in 
a narrow (both in radius and in latitude) boundary layer of the disk which is formed near the 
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contact boundary between the disk and the star. The boundary layer is significant not only due 
to contribution to the luminosity of the compact object. It is important since the boundary-layer 
area is small compared to that of the disk, the layer's emission is harder and must be observed in 
a different spectral range relative to the disk emission. 

Accretion onto a neutron star differs in many ways from accretion onto a white dwarf. 
First, the percentage of energy which is released in the boundary layer can exceed appreciably 
50% because of the general-relativity effects (Sunyaev and Shakura 1986). The boundary-layer 
luminosity Lsl is greater than the disk luminosity even for rapidly rotating neutron stars, 
but the ratio Lsl/ La decreases with the increasing angular velocity of the star (SibgatuUin and 
Sunyaev 1998). Second, at luminosities L > Le^d/lOO ^ 10^^ erg s~^, the local flux of radiative 
energy which is emitted from the boundary layer is of the order of the Eddington flux, i.e., the 
force of radiation pressure on an electron is of the order of the gravitational force acted on a 
proton. There are also other dissimilarities. 

Most of the authors who considered disk accretion onto stars (Shakura and Sunyaev 1973; 
Lynden-Bell and Pringle 1974; Pringle and Savonije 1979; Papaloizou and Stanley 1986; Popham 
et al. 1993), white dwarfs (Tylenda 1981; Meyer and Meyer-Hofmeister 1989; Popham and 
Narayan 1995), and neutron stars (Shakura and Sunyaev 1988; Bisnovatyi-Kogan 1994) suggest 
that the accrcting-plasma velocity decreases from the Keplerian velocity to the stellar rotation 
velocity through the turbulent friction between the differentially rotating layers of accreting 
matter within a thin boundary layer of the disk (sec Fig. 1). In this case, duj/dr > in the 
boundary layer, the boundary-layer meridional extent is ROu = R^2Thi/mp/v!^, and its height is 
hbi ~ R{2Tbi/mp)/{v^)'^, where R is the stellar radius, T^i is the boundary-layer temperature given 
in energy units, rup is the proton mass, and is the Keplerian velocity at the radius R (Shakura 
and Sunyaev 1988). The meaning of the geometrical quantities hu and 6bi is clear from Fig. 1. We 
see that in this approach, hu is determined by the scale height in a gravity held with a free fall 
acceleration go = ^M/R? = {vj^)^/R. The boundary-layer meridional extent, R9bi: is given by the 
same formula that is used to calculate the disk thickness, i.e., only the tangential component of 
the gravity force is taken into account. The thickness, h^i, and the meridional extent turn out to 
be independent of latitude and accretion rate, respectively. 

In this paper, we take a completely different approach to the problem. The deceleration 
of rotation is considered simultaneously with the spread of accreted matter (see Figs. 2-4). 
Discarding popular beliefs, we assume that the deceleration of rotation due to the gripping on 
(hooking, catching on) the slowly moving dense matter beneath the underlying surface, S, is the 
principal friction mechanism. This gripping is due to the turbulent viscosity. This problem has 
much in common with the well-studied (both theoretically and experimentally) problem of the 
deceleration of a subsonic or supersonic flow above a solid surface ( see Subsec. 2. 5 for a detailed 
discussion of the turbulence description accepted in this paper). 

Here, we ignore the magnetic field, H, of the neutron star. As we show in Subsec. 4.8, this 
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Fig. 1. — (a) Shape of the boundary layer in the standard approach. The notation is as follows: e is 
the equatorial plane, D is the disk, and S is the stellar surface, (b) Angular velocity u = Vipjr versus 
radius, where v^p is the rotation velocity; 1 - the Keplerian dependence: u oc r~^/^; the angular 
velocity U)n in the neck is at a maximum and approximately equal to the Keplerian velocity, fn = 
c^„/27r « v^/2ttR = IM{M /Mr^"'!'^ / rI''^ kHz, Rq = i?/(10 km); inside the star from the center 
up to its surface r < R, the rotation is rigid: lo = cos- 
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Fig. 2. — Rotation of the matter in the disk and on the stellar surface in the proposed model; S is 
the stellar surface, p is the pole, and e is the equator. 

assumption is valid for H < 2Ax 10^ {Lsl / LgddT'^'^ gauss, where Lsl and Lf,dd are the spread-layer 
luminosity and the Eddington critical luminosity, respectively. In addition, when deriving the 
equations of motion of the accreting flow over the stellar surface, we disregard its rotation and 
general-relativity effects. An appreciable fraction of the X-ray bursters are neutron stars with a 
rotation frequency fns = uJnsf^T^ of the order of 300 Hz see, e.g., a review by Van der Klis 1998). 
At the same time, the Keplerian rotation frequency fk = cok/^n at the disk-star boundary in the 
problem under consideration is about 2000 Hz. Thus, the effects of stellar rotation, which are of 
the order of the ratio (fns/fk)^, are small, and the computations which are performed below also 
yield an acceptable first approximation to the problem of accretion onto a rotating neutron star. 



1.1. Comparison of the Standard and Proposed Approaches 

The most important points (I-IV) that pertain to our comparison of the approaches are listed 
below 

(I) Let's consider the relationship between the radial (inflow to the star, Vr < 0) and 
meridional {vg) velocities. The meridional velocity is commonly ignored {\vr\ ^ vg). The 
one-dimensional (1-D) approach to the accreting disk also extends to the boundary layer. In any 



-5- 




Fig. 3. — Proposed model: (a) the star is the sphere S, the one-dimensional spread layer is lDi,i, the 
two-dimensional transition between the disk and the spread layer is 2D, and the one-dimensional 

disk is IDd; (b) the linear rotation velocity versus latitude. The rotation is "localized" in 
the latitude belt < < 0^ - for a nonrotating star, the matter outside this belt is essentially 
nonrotating; (c) local radiative flux q versus latitude. The arrows a and b "fence" the ring belt of 
enhanced brightness. In Fig. 3a, this belt lies between the latitudes a and b. 

case, vg is generally omitted from the equations. The boundary layer turns out to be an extension 
of the disk in the sense that the ID description of the disk (IDd) also applies to the boundary 
layer. The disk thickness h passes through a minimum near the neck (see Fig. la). After the neck, 
the flow broadens as it approaches the stellar surface and 'digs itself into the star. 

Meanwhile, a flow with an approximately ID disk, with a 2D region of a sharp turn of the 
flow near the neck, and with an approximately ID (hydraulic) spread of the boundary layer 
over the neutron-star surface seems more natural (see Figs. 1, 3, 4 and Subsec. 2.1). In this 
case, vq 3> \vr \ in the boundary layer. We specify inflow conditions at the interface between the 
boundary layer and the 2D region when calculating the former, because we do not consider here 
an exact 2D problem of the flow in the intermediate zone. In this formulation, it would be more 
correct to talk not about the boundary layer of the disk but about the spread layer. The dynamics 
of the latter bears no relation to the disk. Here, the problem of the deceleration and spread of the 
rapidly rotating matter arises. 

(II) We consider the deceleration of rotation and the structure of isolines of rotation velocity 
Vip. In the standard and proposed approaches, we have Vy, = v^{r) and = v^p{0) in the boundary 
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Fig. 4. — Spread of the rotating plasma from the disk, D, over the neutron-star surface, S. Here, 
/ is the intermediate zone near the disk neck, {) < 6 < 9* \s the hot beh, and 9 > 9* is the cold 
part of the spread layer. The rotation velocity (filled circle) is directed along the normal to the 
plane of the figure. The slowly circulating (in ip and 9) dense underlying layers of matter beneath 
the spread layer are indicated by the dashes. 

layer, respectively. The lines of constant (isolines) are perpendicular to the equatorial plane in 
the standard model and perpendicular to the surface S in the proposed model. In both cases, the 
rotation is Keplerian in the disk at r > r„, where r„ is the neck radius. 

The Vip isolines arc indicated in Fig. 3a by the dots. They are roughly perpendicular to the 
surface S inside the boundary layer. The isolines make a sharp turn near S and then run in a 
dense bundle beside S for a star at rest. This implies that the rotation velocity v,p = ujr abruptly 
decreases near the interface between the spread layer and the underlying surface. The deceleration 
near the surface S is discussed in Subsec. 2.5. 

(Ill) In the standard models, the boundary layer and main energy release concentrate near 
the equator (see Fig. 1). In the proposed model, the width of the rotating belt is variable. It 
increases with increasing accretion rate, M. As M — Mpoie, the rotation extends to the entire 
stellar surface: 0* — > 90°. Calculations show that the spread-layer luminosity which corresponds to 
this accretion rate is Lpgi^ ~ O-QL^dd, where 

Ledd = ^TTR'qo = 1.26 X lO^^^, erg/s, 

qo = = St5oc = 1.00 X 10^5 erg/s cm^ (1.1) 

or Rq 

where go = ^M/B? = {v^y/R, is the Keplerian velocity at the stellar surface, 
St = rUp/aT = 2.5 g/cm^ defines a column of matter with a unit optical depth for 
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Thomson scattering, uip is the proton mass (the plasma is assumed to be hydrogenic), 
(Tr = (87r/3)(e^/meC^)^ is the Thomson scattering cross section, and -Re = R/{W km). Below, 
we talk only about the spread-layer luminosity at the stellar surface normalized to the critical 
Eddington luminosity. Clearly, the accretion disk makes a comparable contribution to the system's 
total luminosity. 

(IV) We consider the latitude dependence of the intensity of radiation q from the belt. In 
the standart approach, the function q{6) has a maximum on the equator. In the approach under 
discussion, q{9) turns out to have a minimum at ^ = 0°. Moreover, q{6 = 0) 0. The distribution 
of q (see Fig. 3c) resembles the letter M, if the plot is complemented with the half which is 
symmetric around the equator. 

1.2. Characteristic Features of the Model 

The problem of the spread of rapidly rotating accreting matter over a neutron-star surface 
has several important features. 

(I) The centrifugal force which acts on the accreting matter in the vicinity of the equator 
essentially offsets the gravitation, because the accreting matter in this vicinity has a velocity that 
is approximately equal to the Keplerian velocity. As one recedes toward the pole, the centrifugal 
force decreases due to the deceleration of rotation. 

(II) In the radiating belt, the radiation pressure works against the difference between the 
gravitational component and the normal component of the centrifugal force. Because of the rapid 
rotation, the local Eddington luminosity at which the radiation-pressure force equals the pressing 
force is appreciably lower than the luminosity qo calculated from the standard formula (1.1). 

(III) The dissipation of rotational kinetic energy causes a strong energy release near the lower 
boundary of the spread layer. As a result, powerful radiation, which diffuses through the spreading 
plasma and forms the observed spectrum, is produced in this sublayer. The normal component of 

the gradient in the radiation pressure (or the normal component of the radiation-pressure force) 
is opposite the direction of the gravitational force. Its direction coincides with the direction of 
the normal component of the centrifugal force. Within the layer, the normal component of the 
radiation-pressure gradient (together with the centrifugal force) decreases the free fall acceleration, 
decreases the density of the matter in the layer, and increases its thickness. An important point is 
that the increase in the height of the column (in the spread-layer thickness) and the drop in the 
plasma density cause the efficiency of the turbulent friction of the layer against the stellar surface 
to decrease. As a result, the deceleration of Keplerian rotation occurs on larger meridional scales 
than that for a small radiation-pressure force. 

(IV) The numerical solution of the equations of radiation hydrodynamics, which is presented 
below, gives the following picture. The local radiative flux, q, at the equator is much lower than 
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the Eddington critical value (1.1). This is attributable to the action of the centrifugal force and to 

the fact that the rotation velocity near the equator is close to the Keplerian velocity. A ring 
zone in which the energy release increases due to the decrease in the centrifugal force is formed 
at a distance from the equator along the meridian. Up to 70% of the layer luminosity is released 
in two such latitude rings (see Figs. 3a and 3c). In this case, the local flux q never reaches the 
Eddington limit (1.1), because the finite rotation velocities Vip are required for the energy release 
through friction against the underlying layers, while the allowance for the centrifugal force turns 
out to be important. The flux q at maximum depends on accretion rate and accounts for 70 — 97% 
of the limit (1.1). 

(V) The spreading plasma within the bright belt is radiation-dominated {pr » Ppi). The 
high saturation of the radiating belt with radiation causes the speed of sound, Cg, in the layer, 
which proves to be much larger than the speed of sound in plasma without radiation, to increase 
significantly. As a result, the flow rotation velocity, v^p, in the layer is moderately supersonic: 
Ma,^ = Vip/cg ~ 5. The meridional spread is subsonic, which is important in choosing a friction 
model. If the Mach numbers, Ma^, were calculated from the plasma speed of sound alone, then 
the flux would be hypersonic with Ma,^ ~ 100; moreover, the flux in would be transonic. 

(VI) At 6 > 9i„ the velocity, u^, drops virtually to zero, and the layer contracts to a thin, 
cold weakly radiating film ('dark' layer), which slowly moves toward the poles. Having lost the 
excess (relative to the star) angular velocity in the radiating belt, the matter becomes cold and 
dense. The intensity of the radiation from its surface decreases by several orders of magnitude. It 
moves relative to the star with a velocity that is considerably lower than the speed of sound. Of 
course, this matter is not accumulated near the poles. There is a slow flow in which it spreads 
under the hot layers over the entire stellar surface. This matter forms the surface against which 
the radiation-dominated layer rubs. 

(VII) If the local radiative flux were the limiting Eddington flux (1.1), the width of the 
radiating belt (let us denote this width by 9^ to distinguish it from 9i,) could be estimated from 
the energy balance 



where Lsl is the spread-layer luminosity. As will be seen from the analysis in the main body of 
this study, (1.1)' underestimates appreciably 9i,. This is attributable to the action of the radial 
component of the centrifugal acceleration. As follows from the calculations which are given below, 
the approximate balance is 



qeffi9) = ^Tgeffie)c, geff{9)=goGeff{9), Geff = l-U''-W^^l-W\ (1.2)' 



where U = vq/v^, W = v^/v^^ is the Keplerian velocity. In (1.2)', we assumed that the 
centrifugal force is almost completely determined by the rotation velocity, v^p^ because the 



(1/2) M {v^S = go 47ri?2 sin 9^ 9^ = &vcsm{LsL/ L^dd) 



(1.1)' 




(1.2) 
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centrifugal force related to the spread velocity is small: i^l/R <^ go- The latitude dependence of 
rotation velocity v^{6) is required to calculate Oi, using (1.2). The approximate validity of relations 
(1.2) and (1.2)' is attributable to the fact that the Eddington radiative flux q^. in the radiating belt 
offsets the difference Qeff between the weight and the radial component of the centrifugal force 
with a high precision [see below Subsec. 1.3 (I)]. The appearance of rings of enhanced luminosity 
follows precisely from this. Indeed, it follows from (1.2)' that q{0) increases with 0, because Vtf,{d) 
decreases with latitude. 

(VIII) We carried out this study in an effort to explain the observed spectra of bursters in a 
wide range of their luminosities. Several variable sources that are X-ray bursters exhibit variations 
in the X-ray luminosity by hundreds of times (see Campana et al. 1998a, 1998b, Gilfanov et al. 
1998). The theory, which we construct here, claims to describe some of the spectral features in the 
light curve of a neutron star as the luminosity changes by two orders of magnitude. 

(IX) The atmosphere under consideration is characterized by the delicate balance of three 
forces: gravitation, centrifugal force, and radiation-pressure force. Under these conditions (just 
as in the atmospheres of massive hot stars), a strong wind which flows away from the zone of 
enhanced brightness is inevitably formed. As a first approximation, we assume below that this 
wind affects only slightly the transport characteristics of the spread layer and its optical depth. 

(X) It follows from the qualitative estimates and the detailed calculations which are presented 
below (see Subsecs. 4.3 and 4.5) that the derived solutions in the radiating belt lie near the 
limit of validity of the shallow-water approximation. The ratios of the effective layer thickness, 
hf^ff = [ I din p/dr I ]~^, in this region (the derivative is calculated near the base of the spread 
layer) to the stellar radius, i?, are typically ~ 0.1. 

Nevertheless, the hydraulic approximation (Schlichting 1965; Chugaev 1982) describes 
satisfactorily the experimental situation at the values of /ig/j/i?. The geometrical characteristics 
of the spread layer are shown in Fig. 5. The spread layer is divided into two distinctly different 
parts at 0^, : the radiating belt (r.b.) and the dark layer (d.L). Their thicknesses, /i^b and hai, differ 
markedly. The thickness h^i is of the order of the scale height in a gravity field with go = jM/R^. 
The slope tan/3 is dheff/RdO. This slope changes only slightly inside the radiating belt (see Fig. 
5a). The angle (3 is very small in the dark layer. The angle [3 in the radiating belt increases 
as the accretion rate, M, decreases. Therefore flows with low M arc described worse by the \D 
approximation. In Figs. 5a and 5b, profiles 1-4 refer to the different values of M. Profiles 1-4 from 
Fig. 5a are repeated in Fig. 5b on a logarithmic scale. This is necessary because the dark-layer 
thickness, h^i, is so small that it cannot be compared with the radiating-belt thickness h^b on a 
linear scale. In addition, the flow separation into two characteristic regions is more clearly seen in 
Fig. 5b. The latitude 6^^ corresponds to the thin transition zone of a rapid change in thickness. In 
turn, the linear coordinates in Fig. 5a are required to show the shape of the radiating belt. 

(XI) The flows in the disk and in the spread layer are closely related. The accreting flow 
is furnished through the disk, while irradiation of the disk from the stellar surface affects its 
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Fig. 5. — Profiles 1, 2, 3 and 4 of the spread-layer (S.L.) effective tfiickness heff{6) for various 
accretion rates AI. The relative luminosities, LsL/Lgdd, that correspond to labels 1, 2, 3 and 4 are 
0.01, 0.04, 0.2 and 0.8, respectively. The stellar mass and radius are 1.4M0 and 12 km. (a) Linear 
scales, the shape of the radiating belt is seen; (b) logarithmic scales, the dark layer and the spread 
layer are seen better than in the upper panel (especially at small M; the scale R6 covers the entire 
distance from the equator to the pole). 

characteristics. This issue is analyzed in Subsec. 2.1, where the disk thickness in the zone of its 
contact with the spread layer is estimated. 

(XII) In a radiation-dominated layer, the hydrodynamic meridional transfer of radiative 
enthalpy toward the poles plays a major role. This enthalpy per unit mass is of the order of the 
gravitational energy ^M/R. 

1.3. Main Results 

We derive the system of equations of radiation hydrodynamics that describes the deceleration 
of rotation of the accreting matter and its meridional spread. The derivation is based on the 
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assumption that Thomson scattering contributes mainly to the opacity and that the energy release 
concentrates in a relatively thin sublayer at the base of the spread flow. The derived equations 
differ markedly from the standard shallow- water equations. A qualitative analysis of the derived 
system of equations and their solutions (Sec. 3) allows us, first, to develop a procedure for 
the effective numerical integration of the equations and, second, to make a clear interpretation 
of the calculations. We carried out a numerical analysis of the system. The problems of the 
numerical analysis, without which an independent verification of the data would be impossible, 
are outlined in Sees. 3 and 4. It is important to set out the procedure for the solution, because 
if this procedure is clear the results can be reproduced. We thus would like to draw the reader's 
attention to the following subsections: (I) smallness of the coefficients, rigidity of the system, steep 
and gentle segments (Subsec. 3.5); (II) multiplicity of the solutions, selection criteria (Subsecs. 
3.6 and 4.1-4.4); (III) circulation of the cooled accreting matter which lost the rotational velocity 
component, its spread and settling (Subsecs. 4.2-4.4); and (IV) modification of the solutions as 
the friction coefficient and the accretion rate are varied (Subsecs. 4.5 and 4.7). 

Below, we dwell on the most important results. 

(I) Rings of enhanced brightness. Typical latitude profiles of the flux q are shown in 
Figs. 6-8. The radiating belt broadens as the accretion rate increases. The latitude Omax at which 
the bright rings lie increases with increasing M and Lsh/Ledd- The bright rings are not narrow; 
their width is /S.9 ~ At a high accretion rate M ~ Mpde [see Subsec. 1.1 (III)], the entire 
stellar surface, except the small equatorial and polar zones, emits almost uniformly. In this case, 
the bright ring spreads out. Profile 4 in Fig. 7 refers to this situation. 
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Fig. 6. — Local radiative flux q{6) (solid lines) and effective free fall acceleration (including 
the centrifugal acceleration) g(,ff{9) (dashed lines) versus latitude. An important point is that 
q/qo ~ Qeff/do ill the radiating belt [see (1.2) and (1.2)']. Labels 1 and 3 refer to LsL/L^dd = 0-01 
and 0.2, respectively. 

The function q{9), which gives the radiative flux, has the maximum at Omax with which the 
ring of enhanced brightness is associated. This ring is marked by arrows a and b in Fig. 3c. The 
ring bounded by latitudes a and h which correspond to these arrows is shown in Fig. 3a. To the 
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Fig. 7. — Changes in the profiles of the local radiative flux q (solid lines) with M. The radiating 
belt broadens with increasing M. The boundary 9^, lies at the right edge of the belt in the zone of 
an abrupt decrease in q{9). The Qeff/go profile at Lsl/L^cLcL = 0.8 is indicated by the dashed line. 
The scale RB terminates at the pole. Labels 2, 3, and 4 refer to Lsh/L^dd = 0.04,0.2, and 0.8, 
respectively. 




M, km 

Fig. 8. — q{6) profiles in a wide range of M. Labels 1, 2, 3, and 4 refer to Lsh/Ledd = 0.01, 0.04, 0.2, 
and 0.8, respectively. 

right of Omaxi the rotation velocity vanishes (see Figs. 3b and 3c). For this reason, q{0) rapidly 
decreases to the right of Omax- The increase in q{9) to the left of 9max (see Fig. 3c) with increasing 

9 is caused by the decrease in v^p (see Fig. 3b). Indeed, the decrease in leads to an increase in 
the acceleration g^ff ~ ^/o (1 — W^), where W = v^/v^, and, consequently, to an increase in q, 
q oc g^ff [see formulas (1.2) and (1.2)']. 

(II) Surface density, optical depth, and spectral hardness. In the steady-state case 
in the absence of sources and sinks of matter in the spread layer, the flux of mass is conserved 
when spread along the meridian. We ignore the settling of matter to the bottom within this layer 
(Subsecs. 4.3 and 4.4). The meridional velocity is thus oc I/S5, where S5 = / pdr is the surface 
density, the integral is taken over the spread-layer thickness. Accordingly, S5 reaches a minimum 
at the maximum of vq. These profiles are illustrated in Fig. 9. The spread layer consists of the 
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Fig. 9. — Profiles of the meridional velocity vg and of the spread-layer surface density for 
Lsh/Ledd = 0-2. At Oi,, the flow abruptly slows down, and its density increases. In the upper 
corner, the Keplerian velocity is shown for M = IAMq and R = 12 km. 

radiating belt and the dark layer which border at 9^,. The maximum of vo and the minimum of E5 
lie near the edge of the radiating belt. In the dark layer, E5 abruptly increases, while vg abruptly 
decreases. This is clearly seen in Figs. 9 and 10. The dark layer terminates at the sonic point. 
This important problem is analyzed in Subsecs. 4.1-4.4. 

Figure 10 shows how the spread-layer surface density S5 depends on accretion rate. The 
spread-layer optical depth tt = S^/Sy for Thomson scattering rapidly decreases with decreasing 
M. At Lsl/ Ledd ~ 10^^, we have tt ~ 3. No blackbody radiation can be produced in the medium 
at such a small optical depth; the protons and electrons have temperatures much higher than 1 
keV. The low-frequency photons emitted by the dense underlying layers are comptonized by hot 
electrons of the spread layer and form the experimentally observed hard power-law tails in the 
spectra (Barret ct al.l992, Campana ct al. 1998a, 1998b, Gilfanov et al. 1998). In turn, some 
of the hard photons penetrate deep into the underlying layer and heat it up [see Sunyaev and 
Titarchuk (1980) for a discussion of this problem]. 

At a high spread-layer luminosity Lgi > O.lLe^^rf, we have a different picture. In this case, the 
spectrum is much softer and exhibits no pronounced hard tail. The optical depth for Thomson 
scattering is great (tt ~ 1000). The layer optical depth for effective absorption ^Jtt t//(z/, T^) 
(Shakura 1972; Felten and Rees 1972; Shakura and Sunyaev 1973; Illarionov and Sunyaev 1972) 
turns out to be also large. Deep in the layer, the most important physical process is saturated 
comptonization. The parameter y = (Te/meC^) Tj^, which characterizes comptonization in a layer 
with an optical depth tt, is large even at tt > \/m(,c^ /T^ ~ 20. At Tg ~ T ~ t^^'^ (see Subsec. 2.2) 
and Tt ~ 10^, we have y ~ {Ts/nieC^) t^^^ ~ 10^ deep in the layer. The dimensionless frequency, 
near which the rate of photon absorption due to the free-free processes equals the rate of photon 
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Fig. 10. — Shape of the T,s{0) profile at different spread-layer luminosities. At a fixed M, the 
surface density rises with increasing 9 at 6^, after which an optically thick weakly radiating layer 
begins. Labels 1, 2, 3, and 4 refer to LsL/Lgdd = 0.01,0.04,0.2, and 0.8, respectively. 

removal upward along the frequency axis via comptonization, is xq = hvo/Tf, = 3 x 10^ nj"^ jxt^'^ . 
This frequency, deep in the layer, is high. Here, the electron density is in cm~^, and the electron 
temperature is in Kelvins. As can be seen from the plots in lUarionov and Sunyaev (1975), 
comptonization produces the Wien spectrum deep in the layer at y > 10, while at y > 10^ 
the spectrum is similar to a blackbody due to the comptonization of low-frequency photons. 
Consequently, the combined effect of bremsstrahlung and comptonization produces a blackbody 
spectrum deep in the layer. The emission from the dense underlying layers is yet another source 
of soft photons for comptonization. Nevertheless, the emergent spectrum differs markedly from a 
blackbody spectrum, because at each frequency we see a portion of the spectrum that is produced 
at its own depth [see lUarionov and Sunyaev (1972) for a discussion of this problem]. 

Integrated spectrum of the spread layer. The spread-layer spectrum is given by the 

integral of the emergent spectra at each point of the layer over its surface. Here, it is important 
that the spectrum and the local radiative flux depend on latitude. The integrated spectrum is not 
presented for the following two reasons. 

(a) A considerable fraction of the spread-layer radiation falls on the disk surface (see Fig. 
2). This fraction is partially reflected from the disk and partially reprocessed by it into softer 
radiation. The intensity of the disk-reprocessed radiation depends on the binary's inclination (see 
Lapidus and Sunyaev 1985). Allowance for these effects requires a special study. 

(b) Because of the shielding by the disk, the observer cannot see the lower half of the neutron 
star. Part of its surface is located in the shadow region and does not contribute to the direct rays 
viewed by the observer (here, we disregard the general-relativity effects). In Fig. 11, the flux of 
"direct" photons from the stellar surface is plotted against the cosine of the binary's inclination 
angle, which is measured from the star's polar axis (see Figs. 1-4). In these calculations, we 
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assumed the indicatrix of emission per unit area to be (^(/x) = 1 + 2.06 /x (/x = cos?); this is the 
first approximation for the angular dependence of the intensity of the emergent radiation from 
a scattering atmosphere (Sobolev 1949; Chandrasekhar 1960). In the calculations, we took into 
account the profiles q{6) (see Figs. 6-8). 




0.2 0.4 0.6 0.8 1 

fj,= COS i 



Fig. 11. — Functions /(/x)//(0) [/(erg s~^ stcradian""^, /x = cosi] which give the intensity of the 
radiation from the visible neutron-star surface. The binary's inclination, i, is measured from the 
polar axis. Labels 1, 2, 3, and 4 refer to LsL/L(.,]d = 0.01,0.04,0.2, and 0.8, respectively. The 
functions are normalized to the intensity of the radiation from the equator /(O). 

Clearly, the presence of a " shadow" and the reflection by the disk must lead to an appreciable 
polarization of the emission from an accreting neutron star. At low luminosities, at which the 
radiating belt is bound by the equatorial zone, the flux toward an observer in the equatorial plane 
{ji = 0) exceeds the flux toward the polar axis (/x = 1). Curves 1-3 in Fig. 11 correspond to these 
luminosities. At high luminosities, the emission extends to the entire stellar surface (curve 4 in 
Fig. 7) . In this case, the flux toward the polar axis is greater than the flux toward the equatorial 
plane (curve 4 in Fig. 11). Indeed, the ratio I(l)//(0) = 2 for a uniform emission of the stellar 
surface, because the disk occults the lower stellar hemisphere. 

(Ill) Deceleration efficiency. The meridional distributions of rotation velocity are 
shown in Fig. 12. We performed the calculations for 5 = 10^^, where 5 = 1 — Wq = {v^ — v^)/vj^, 
where is the initial rotation velocity in the spread layer when passing from the disk to the 
layer, and is the Keplerian velocity (see Subsecs. 3.4 and 4.6). Similar results are also obtained 
for different values of (5 ^ 1. As M increases, the deceleration of rotation requires an increasingly 
large area of contact between the radiating belt and the stellar surface. The boundary of this belt 
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R6. km 

Fig. 12. — Meridional profile of the rotation velocity. The radiating belt is simultaneously the 
rotating belt (cf. Figs. 6-8). Labels 1, 2, 3, and 4 refer to Lsh/Ledd = 0.01,0.04,0.2, and 0.8, 
respectively. 



An interesting characteristic of the deceleration efficiency is the number of rotations a test 
particle makes around the star before it looses its rotation velocity and reaches the latitude Oi,. 

Table. Dependencies of N^ot and q/Q^ on Lshl^^edd- 



Lsh/Ledd 


0.01 


0.04 


0.2 


0.8 


Re^, km 


0.43 


1.36 


4.5 


16.4 


Nrot 


0.3 


1.1 


5.4 


36 




1.01 


1.05 


1.3 


3.8 



The table shows how the width of the belt 9^, and the number of rotations, 

Nrot = A^/2n, A^= ujdt^ -, 

Jeo Jo V0 COS 6* 

of its rotating plasma depend on the spread- layer relative luminosity LgL/L^dd. 

As we see from the table, at a low spread-layer luminosity (O.OlLedd), a test particle does not 
make even one rotation around the star, rising by 430 m along the meridian and co vering 23 km 
in latitude. Note that the meridional boundary-layer extent (see Fig. 1) ROu = RJITyilm^jv^ 
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is a mere 70 m in the standard approach (at T^i = 3 keV). At a high luminosity, the deceleration 
efficiency is much lower: at Lsl = O.SLgcid, test particles make 36 rotations around the star before 
they slow down and reach 6^. 

(IV) Advection of radiative energy. The spreading matter in an optically thick layer 
transports radiative energy along the meridian from the equatorial region into the region of the 
bright rings. The energy transfer is described by the energy equation (see Subsecs. 2.7 and 3.2). 
Our calculations (see Fig. 13) show that the energy flux q emitted per unit area of the spread 
layer in the equatorial region is much smaller than the local energy release through the turbulent 
friction in a column Q'^ = J etdy, where the integral is taken over the spread-layer thickness 
[Subsec. 2.5, (2.17)]. On the other hand, the flux q at higher latitudes for larger luminosities 
exceeds appreciably [see the table, which gives the dependence of maximum iq/Q^)max in the 
spread layer on Lsl/ Ledd\- This is only to be expected, because the spread-layer optical depth 
for Thomson scattering is large at high luminosities. As a result, the characteristic diffusion 
time of the photons from the base of the hot layer turns out to be of the order of the time of 
hydrodynamic motion through the radiating belt, where ^ {h/c)TT and ~ R9i,/vQ. The 
advection of radiative energy favors the formation of bright latitude rings at the stellar surface. 

8, grad 

10 20 30 



1 

+ 

\ 

0.5 


12 3 4 

Re, km 

Fig. 13. — Spread dynamics is determined in many ways by the hydrodynamic transfer of radiative 
energy. The ratio of the energy flux, emitted per unit area of the radiating layer to the frictional 
energy release per unit area of the contact surface between the spread layer and the star. The 
energy is transferred from zone A into zone B through meridional advection. This calculation is 
for Lsh/Ledd = 0.2. 



1.4. Problems That Require Further Studies 

(I) Possible Instability of the Spread Layer. The energy-producing sublayer props up 
the spread flow from below, and the radiation pressure dominates. This may be the cause of 
convective instability with floating 'photon bubbles'. Clearly, this requires that the spread layer 




- 18 - 



be optically thick. The fluctuations which are produced by the bubbles may be responsible for 
the observed low-frequency noise and even for quasi-periodic oscillations in accreting neutron 
stars. Although, on the other hand, the rotation velocity, v^, inside the spread layer outside the 
energy-producing sublayer, along with the normal component of the centrifugal force, slightly 
increases with height. In addition, the tai]|] of the volume density of dissipative energy production 
it erg s~^ cm~^ extends from the energy-producing sublayer into the spread layer. As a result, 
the flux q and, hence, the Eddington force slightly increase with distance from the surface. These 
two effects play a stabilizing role with respect to the convective instability. The instability is also 
stabilized by turbulent viscosity and turbulent diffusion, which are attributable to the presence of 
a wall (the stellar surface from below) and of shear turbulence. This turbulence causes, flrst, the 
mixing of density fluctuations (diffusive mixing of photon bubbles) and, second, the momentum 
exchange between the bubbles and the ambient medium (viscous deceleration). The problem of 
the vertical stability of hydrostatic quasi-equilibrium, indubitably, requires further study. 

(II) Gap between last stable orbit and surface. Here, we do not consider the case where 
the stellar radius R is smaller than the radius of the marginally stable Keplerian orbit Rc = 3Rg. 
In this case, part of the gravitational energy is released in the zone of contact between the disk 
and the star due to an appreciable radial velocity (it increases with the increasing R^ — R) of 
the accreting matter on the stellar surface. The energy release at the equator causes the layer to 
swell in this zone because of the radiation pressure, and the ID approximations appear to become 
inapplicable. 

(III) Rapid rotation. The case of rapid rotation of a star, where its angular velocity, ojs, 
accounts for an appreciable fraction of the Keplerian angular velocity, uj^i requires a special 
analysis. 



1.5. Structure of the Paper 

Here, we derive and numerically solve the system of equations of motion of the matter over 
a neutron-star surface in the ID approximation. The vertical (along r) and horizontal (along 
the meridian 6) scales are separated in this approximation. The gradients in r are assumed to 
be large compared to the gradients in the polar angle 9. For this reason, the vertical profile is 
constructed from the conditions of hydrostatic quasi-equilibrium. The profile is given by simple 
analytic expressions. The averaging of the quantities over the vertical structure with which the 
derivation begins allows us to obtain a system of equations for these quantities that contains only 
gradients in 9. These are the equations of radiation-dominated shallow water that describe the 
dynamics of the high-velocity layer in which the newly accreted matter fiows over the surface of 
the slowly moving high-density underlying layers (Sec. 2). This completes the part of the paper 



The frictional heat release concentrates mainly in the energy-producing sublayer; see Subsec. 2.5. 
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in which we formulate the problem. The part in which we present the results is devoted to the 

steady-state solutions of our model. In this part, we describe the physics of adaptation: effects of 
the friction model, Eddington swelling of the boundary layer, regulation of the deceleration by 
the swelling, etc. (Sees. 3 and 4). We study the mathematical properties of radiation-dominated 
shallow water: the structure of the solutions, the limiting set of steady-state solutions, its 
one-parameter structure, the critical solution, etc. (Sec. 3). We analyze the dependences of the 
boundary-layer characteristics on M and on model parameters (Sees. 3 and 4). Note that it 
follows from the solutions that the bulk of the radiating belt is roughly in hydrostatic equilibrium 
along the meridian, as was previously suggested by Shakura and Sunyaev (1988). In this case, 
the "pressing" force to the equator is the tangential component of the centrifugal force, which 
counteracts the meridional component of the pressure gradient. 

The interaction of the spreading flow with the underlying matter is described in Subsec. 2.1. 
Here, we also estimate the effects of disk heating by the spread-layer radiation and determine the 
minimum possible disk thickness near the neck. The neck forms a kind of base of the disk at the 
location of its interaction with the spread layer. In Subsec. 2.2, we determine the structure of a 
radiation-dominated atmosphere in which the scattering by free electrons is a major contributor to 
the opacity. The force and energy characteristics of this atmosphere which arc required to calculate 
its energy and transport properties are derived in Subsec. 2.3. A simple formula for radiative 
losses in the spread layer is written out in Subsec. 2.4. Subsec. 2.5 deals with the friction model. 
In Subsecs. 2.6 and 2.7, we derive the equations of "radiation" shallow water. Transformations 
of this system of equations are made in Subsecs. 3.1-3.3. In Subsecs. 3.4 and 3.5, we study the 
phase space of the system and introduce an important surface - the "levitation" surface. The 
dependence of the spread velocity on the Mach number is analyzed in Subsec. 3.6. The solutions 
of the derived system of equations are analytically and numerically analyzed in Subsecs. 4.1 and 
4.2. Subsec. 4.3 investigates the critical regime in which the spread velocity reaches the speed 
of sound for the spread of matter at the edge of the radiating belt. The dynamics of the matter 
spread outside the radiating belt is analyzed in Subsec. 4.4. The effect of friction on the results 
is considered in Subsec. 4.5. Such an analysis is necessitated by the fact that, first, the friction 
coefficient is determined only approximately and, second, this coefficient decreases by a factor of 

1.5 to 2 in the supersonic regime (see Subsec. 2.5). We construct the dependence q{9) in Subsec. 

4.6 and study how the spread flow is modified as the accretion rate M is varied in Subsec. 4.7. 
At large M, the "squeezing" of the flow due to the convergence of meridians to the pole becomes 
important. In addition, the spread flow in the radiating belt is strongly subsonic at high M and 
transonic at low M. In Subsec. 4.8, we estimate the limiting magnetic field strength that can 
still be ignored when the deceleration and spread are studied. The results are summarized in the 
conclusion. 



2. HYDRAULIC APPROXIMATION 
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2.1. Spread Scheme 

The spread flow is illustrated in Fig. 4. Here, e is the equatorial plane, D is the disk, / is the 
transition region between the disk and the boundary layer (the region of 2D flow), and S is the 
stellar surface. Let h{0) be the boundary-layer thickness. For the ID approximation to be valid, 
it is necessary that 

/i < i?, \dh/dl\ < 1, (2.1) 

where I = RO is the arc length along the meridian. Note, incidentally, that it follows from 
the calculations (see below) that the angle arctan |d/i/dZ| is constant in order of magnitude at 

The latitudinally and mcridionally moving spread layer draws (grips) the thin underlying 
layers into motion via turbulent friction. The underlying layer is indicated by the dashed line in 
Fig. 4. The distribution inside the boundary layer is roughly uniform in height (see Subsec. 
2.5 for a discussion). Almost all of the drop in takes place near the surface, 5, which is the 
bottom of the boundary-layer, r = rs oi y = ys- The main turbulent energy release occurs near 
the bottom. Since the problem is a steady-state problem, the temperature T is constant at y <ys 
(these layers were heated earlier, and the energy release in them is small), while the vertical 
radiative flux q is approximately constant above ys and equal to zero below. Because of the abrupt 
change in q at the surface S, the plasma parameters also change when this surface is crossed. The 
underlying-layer density p2 exceeds the boundary-layer density in the hot subregion < 9 < 9^,. In 
the region 9 > 9-/,, the density of the underlying layers rapidly increases with depth, because they 
are isothermal. For this reason, the velocity in the underlying layers rapidly decreases with depth. 

Disk heating by spread-layer radiation. Let us find out how the disk thickness near a 
neutron-star surface will change under spread-layer radiation. The angular rotation velocity uj of 
the accreting matter is known to reach a maximum in the transition zone between the disk and the 
spread layer (Shakura and Sunyaev 1973). Accordingly, duj{r = rn)/dr = 0. This point in the flow 
is called the neck (see Figs. 1, 3, and 4). The dissipative production of thermal energy vanishes at 
r„, because the derivative of oj{r) at this point changes sign and because the heat release, which is 
proportional to (du/dr)^, is small in the neck zone. The surface density of viscous energy release in 
the disk (in Newton's approximation) behaves as = (37MM/47rr^)(l — ^jR/r). At r » i?, the 
energy release is a factor of 3 greater than the local gravitational energy release ^MM /Airr^^ 
which is attributable to the transfer of mechanical energy by viscous forces over the disk from 
regions with smaller r to regions with larger r. At the same time, as r ^ R. Therefore, the 

disk half-thickness, h = (3/87r) (ctt M/nip c)(l — ^jR/r), also decreases as one approaches the neck 
(Shakura and Sunyaev 1973). This consideration ignores disk heating by the spread-layer radiation 
and the pressure of external radiative flux. Let us estimate how the heating by external radiation 
affects the disk half-thickness near the neck. Let T^, Tc, and be the disk-surface temperature, 
the temperature in the disk equatorial plane, and the external-radiation temperature, respectively. 
Since the flux from the spread layer accounts for a fraction of the Eddington flux, its temperature 
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isTr = 1 — 1.5 keV. The disk albedo is assumed to be 0.3-0.5. Consequently, the disk surface near 
the neck will be heated to temperatures of the order of 1 keV. Clearly, the surface temperature 
under quasi-steady-state conditions (the radial velocity of the matter near the neck is low) is less 
than or equal to Tc- Using a hydrostatic estimate of the half-thickness, h = Rcs/vj^, and assuming 

that Cg ~ '^Tc/rup, we obtain h = 30Rq^^ ^/Td \J M jM^ m, where is in keV, and R% = R/{10 
km). It thus follows that the disk angular half-thickness, is finite. The conclusion about the 
finite but small disk half-thickness near the neck is crucial for the validity of our approach to the 
problem of accretion onto a neutron star. 

Here, we do not consider the advection solution for the flow in the disk (Narayan and Yi 
1995) for two reasons. First, the neutron star has a solid surface, and it is, therefore, difficult to 

imagine accretion without energy release on the surface. Second, the flux from the spread layer is 
so large that it will cool the plasma flow at a distance of several tens of neutron-star radii through 
comptonization. As a result, an accretion pattern with a relatively thin (geometrically) accretion 
disk is most likely realized near the star. 



2.2. Structure of the Spread Layer (Atmosphere) 

We make the following assumptions: 

(I) The energy is mainly released in a thin sublayer near the boundary between the spreading 
and underlying layers. 

(II) Thomson scattering makes a major contribution to the opacity. 

(III) We define the free fall acceleration g^ff within the spreading layer with allowance for the 
contribution of the centrifugal acceleration as the difference ^g// = ^M/E? — v'^/R — v^/R [see 
(1.1) and (1.1)']. For completeness, we added the centrifugal acceleration due to the meridional 
motion. It is small compared to the acceleration ^M/R^. The rotation velocity, v^, is assumed to 
be constant within the atmosphere (if is fixed). 

Under these assumptions, the atmospheric structure is described by the system of equations 

P = 9eff S, (2.2) 
cTtT de 



3 dS' 



(2.3) 



P = 2^ + |, (2.4) 

e = aT^. (2.5) 

The system of equations of hydrostatics (2.2), radiative heat condiiction (2.3), thermodynamic 
state (2.4), and Stefan-Boltzmann's law (2.5) relates the unknown functions p,e,p, and T of the 
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vertical column density S, which is a Lagrangian coordinate {dTi/dy = p). The y and E axes are 
directed downward. Accordingly, the pressure in (2.2) increases with S. The q is the absolute flux. 
We assume that the layer borders vacuum (p = 0) aty = 0,S = 0. The boundary-layer bottom 
y = ys = h, S = S5 is the surface S (see Fig. 4). The system of equations (2.2)-(2.5) is written 
for an arbitrary fixed latitude 6. Due to the approximate assumption^ of constant v^{9) in the 
boundary layer < y < h, the acceleration g^ff in (2.2) is constant as S increases deep into the 
layer. 

Eq. (2.3) holds for radiative transfer in an optically thick plasma with ^ Tit with 
dominating Thomson scattering. We assume that the turbulent dissipation is localized near the 
surface S. It thus follows that the heat fiux q in the spread layer < y < h, which is transferred 
via radiative heat conductance, is approximately constant. The total pressure consists of the 
plasma and radiation pressures. 

The system of equations (2.2)-(2.5) can be easily solved. The solution that satisfies the 
boundary condition 

p{0) = 0, e(0) = 0, p(0) = 0, r(0)=0, (2.6) 



IS 



V 3g rup g^r ( ac ^V/"^ f3q 
P = 9effJ^, e=-rT, p= [-r^) , ^ = (^- j , (2.7) 

dwr = 90 Gyjr, Gwr = A = Gg// — G,-, Gr = — , = ^ j (2-8) 

go cL,T 

where Ge// = geff/do [see Subsec. 2.2 (III)]. The difference. A, of the gravitational acceleration, 
the component of the centrifugal acceleration normal to S, and the radiation acceleration gr (2.8) 
is particularly important for the subsequent discussion. 

Let us derive the dependence on y. The Lagrangian and Eulerian differentials are related by 
dT, = pdy. Substituting the solutions (2.7) and (2.8) into this relation and integrating it using the 
conditions (2.6), we obtain 

■y =-r^ P = o oin y ' I3 = rn g^acET, (2.9) 



Ppl = Qwr S, Pr=gr^, P = Ppl + Pr = Qeff S, T = y. (2.10) 

o 

These solutions are very simple. The atmosphere under consideration is similar in structure to 
the atmospheres of X-ray bursters during outbursts and to the atmospheres of supermassive stars. 
We describe these solutions in some length, because the calculations of transport characteristics 
are based on them. 



is assumed to depend on latitude alone. By we mean the r-averaged velocity in the boundary layer. 
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2.3. Lateral Force and the Surface Energy Density 

Here we study the dynamics of a layer of shallow water or a fluid film that flows over the 
underlying surface S (see Fig. 4). The thickness, h, of the layer is assumed to be small compared 
to its horizontal extent R9 [see (2.1)]. When the equations with gradients in 9 are derived, the 
layer is broken up into the differential elements R59. The force interaction between adjacent 
elements is effected by the lateral force pdy. Substituting (2.9) and (2.10), we obtain 

/ pdy = ^—(^ , 2.11 

Jo 5 A nip 

where, as above, the subscript S denotes values near the bottom of the spread layer. It is 
important that the force (2.11) contains the large dimensionless factor Ge///A at A ^ 1. 

The layer transports mass, momentum, angular momentum, and energy over the spherical 
surface S. In particular, the advection of radiative energy takes place. The expressions for 
the surface densities of mass, momentum, and angular momentum are obvious: E5, vgT,s, 
R cos9v^ Tis, where vg{9) is the r-averaged meridional velocity in the layer, and Rcos9 is the arm 
relative to the rotation axis, which the polar axis is. Let us calculate the surface energy density 
(erg cm^^) using the distributions (2.9) and (2.10). The total energy of the layer consists of the 
kinetic energy K = {vg + v^) S5/2, the gravitational energy {Eg), and the internal energy (Eint)- 
In turn, the latter consists of the plasma (Epi) and radiative (Er) contributions. The surface 
density of gravitational energy is 



E, = t P9eff {h - y)dy = | ^ (2.12) 
Jo o Z\ nip 



Interestingly, the expressions (2.11) and (2.12) are identical. The energy (2.12) is measured 
from the bottom of the layer y = h. Since the volume density of plasma internal energy is 
3pT/mp = {3/2)ppi, we have Epi = (3/2) Jq Ppidy. Calculating this integral and the integral 



Er= r edy = (24/5) (G,/A) Tg ^s/nip, 
Jo 



we obtain 

12 Ge// + Gr Ts Sg 

- y — A ^^-^^^ 

All expressions (2.11)-(2.13) that appear in the formulas for horizontal interactions in the layer 
increase as the layer swells with A — 0. 



2.4. Radiative Cooling 



The energy losses by the layer surface are determined by the flux which comes from the inside 
and can be readily expressed in terms of the layer surface temperature, T^, on the side of vacuum 
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{q = acT^/4). This is the temperature at an optical depth of the order of unity (S^ ~ Sy). The 
solutions (2.7) and (2.8), in which the integration constant for the equation of radiative heat 
transfer (2.3) was discarded, is valid at large S (S ^ '^t)- Formula (2.7) for T allows the local 
flux, q, to be expressed in terms of the local S and T [q = acT'^ /3tt)- In particular, inside the 
spread layer near its bottom, we obtain 

g = acTj St/SSs. (2.14) 

2.5. Turbulent Deceleration of Rotation and Viscous Energy Release 

Profiles of average and fluctuation velocities. Let us dwell on friction and dissipation 
and analyze the friction of the boundary layer against the star. The turbulence via the friction it 
produced decelerates and heats up the layer. In this case, the rotation and meridional spread of 
the matter in the layer draws the underlying layers into motion through turbulent viscous gripping. 
This issue is discussed in Subsec. 2.1. Since the source of heat concentrates near the boundary S, 
the temperature under the surface S is constant, while the entropy s rapidly decreases with depth. 
Stratification with a large reserve of gravitational stability arises (the Richardson number is great 
in the underlying layers). It is difficult for the high-entropy layer, which moves from above, to 
"grip" the low-entropy underlying layer. This situation roughly corresponds to a wind above a 
smooth surface. 

Turbulence in a gaseous layer that flows with a subsonic velocity above a flxed boundary 
has been well studied and described by Loitsyanskii (1973), Schlichting (1965), Chugaev (1982), 
Landau and Lifshitz (1986), and, recently, by Lin et al. (1997). For low viscosity, the height 
distribution of the mean velocity is described by Prandtl-Karman's universal logarithmic profile 

{v) ^^v^ ln(7.7 Ay/i^), r = pvl (2.15) 

where Ay is measured from the surface S, and r is the tangential stress which does not depend on 

Ay. The parameter is important. It determines the turbulent velocity and pressure pulsations. 
The universal profile is such that is velocity pulsations about the mean in all directions: ~ 

V {{ve - {vem ^ VW) ^ ^| {(v^ - {v^m or 

V* = \J {{v- (v^) 4 - (ve) 4 )2 ) /3, 

where v = {vg,Vr,v^}, and x,z is the wall plane. This implies that the pulsations are isotropic 
at a given distance. Ay, from the wall in a comoving frame of reference which moves with the 
mean velocity {v) of the flow at this distance. The order-of-magnitude amplitude of the pressure 
pulsations about the mean is pv^. An important point is that the velocity -y* and the tangential 
stress r in the region of turbulent mixing do not vary as the distance to the wall Ay is varied. 
Actually, the law (2.15) follows from the condition of the invariance of r. When this law is 
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deduced, the stress is written in Newton's form r = putd{v)/dy, and the mixing-length theory 
is used, in which the turbulent viscosity ft = It v-^ is determined by the scale of the wall vortices 
If = KAy. The coefficient k ~ 0.4 is called Karman's coefficient. The integration constant in 
Newton's friction law {v) = J dyr . . . , which appears in (2.5) under the logarithm, is chosen by 
joining the viscous and turbulent solutions in the transition region between the viscous sublayer 
and the turbulent region. Formula (2.15) holds outside the viscous sublayer, when Ay > 1^, 
v-^lyjv K, 12, where \y is the thickness of the viscous sublayer which is determined by the molecular 
and radiative viscosities. 

The disk-accretion theory commonly uses an expression for the kinematic viscosity ft of the 
form vt = ahd {cs)d, where hd and {cs)d are the disk thickness and the speed of sound in the 
disk, respectively. For the Keplerian dependence of angular velocity tOd on r, this expression 
reduces to a simple formula for the stress r = ap, where p is the pressure in the disk (Shakura 
and Sunyaev 1973). Let us compare the formulas Ut = ahd {cs)d and ft = nAyVi,. We see that 
the wall turbulence differs from the disk turbulence, first, by the form in which the mixing length 
is written {k Ay instead of hd) and, second, by the formula for turbulent velocity pulsations [v* 
instead of a {cs)d\- 

Friction coefficient. Let us calculate the stress r. For this purpose, we extend the 
distribution (2.15), which is weakly dependent on Ay, to the height Ay = h. To obtain an 
estimate, we assume that at this height the mean velocity of flow {v) is equal to the Keplerian 
velocity. The ion viscosity is fj = 2.2 x 10~^^ T'^'^ / ps InA = 1.3 x 10^ cm^ s~^, where, T is in 
degrees, and ps is in g cm"^ (Spitzer 1962). We also assume the following: h = 100 m, S5 = 300 
g cm~^, T = 3 keV, and the Coulomb logarithm In A = 10. This viscosity is small compared to the 
radiative viscosity 

4 a T'^ St Pr 
f r = TT 5— ^ — It c 

15 p 

(Weinberg 1972); here, pr = aT^I(? is the photon "density". It = l/nax is the photon mean 
free path, and n = p/mp. It follows from the numerical calculations that the radiation density pr 
accounts for a small fraction (~ 10~^) of the plasma density ps- For the above values of p and T, 
we obtain = 1.7 x lO''' cm^ s~^. The mixing is governed by the radiative viscosity: ^ z/j. 
Substituting f,. in (2.15) and setting M = 1.4M0, and R = 12 km when calculating the Keplerian 
velocity, we arrive at the equation 4.2 x 10*^ = ln(';;^/220), where the unknown is in cm s^"*^. The 
estimates are given in the approximation of Newton's potential. Solving this equation, we obtain 
3 x 10^ cm s"-*^. In this case, the viscous-layer thickness 1^ ~ 12f /v* 0.7 cm is much smaller 
than the scale height h. The Reynolds number is Re= hv^/u ~ 10^. Let us write 

T = abpv'^. (2.16) 

The coefficient of friction against the bottom is ai, = v^jv^ [cf. (2.15) and (2.16)]. In our example, 
OLh ~ ~ 0.81 X 10~^. The value of changes only slightly as the parameters are varied. 

In Subsec. 3.1, we compare the terms which are associated with a}, (wall turbulence) and a (a 
coefficient that is widely used in the disk-accretion theory). 
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Dissipation and the height distribution of energy release. The volume power density 
(erg cm~^) of the heat source et attributable to turbulent friction is given by 

et{Ay) « (5/2) r/Ay (2.17) 

(Landau and Lifshitz 1986). It follows from (2.17) that the surface density of the total viscous 
energy release is 

Q+ = et{Ay) d{Ay) « et d{Ay) ^ 2.5tv^ In ^ « tv. (2.18) 

Jo Jl^ Iv 

Here, v is the height-averaged velocity in the layer. In the case under consideration, whe re the 
flow moves latitudinally at velocity v^p and meridionally at velocity vq, the total velocity + v'^ 
must be substituted in (2.18). The energy release (2.17) takes place mainly near the bottom of the 
boundary layer. For example, 

[/ etd(Aj/)]/Q+«0.2. 

Jo.lh 

An analysis of the calculations which are given below shows that Vip is several times faster than 
the speed of sound Cg at 9 < 9-), - the flow is moderately supersonic. In this case, the expressions 
(2.15)-(2.18) remain valid (Loitsyanskii 1973; Schlichting 1965). In a supersonic case with the 
Mach number Ma ~ 5, in (2.16) decreases by a factor of 1.5 to 2. 

Above, we described the surface characteristics (Subsec. 2.3), radiative cooling (Subscc. 
2.4), and dissipative heating (Subsec. 2.5). Now, we can proceed to a derivation of the balance 
equations that relates to all these physical processes. 



2.6. Equations of mass, momentum, and angular momentum transfer 

Let us derive the first three equations of the system of four ID equations which represent the 
laws of conservation of mass (equation S), momentum (or angular momentum along the meridian) 
in 9 (equation vq), angular momentum (equation v^), and energy (equation T). We take a film 
annulus with the area dS = 2ttRC R59 on the sphere S, where C = cos 9. Its mass is dS S5. The 
mass flux through the 9 = const conic surface (annulus boundary) is 2TrRCvQTis- Accordingly, the 
equation (E) is 

RC{J:s)i = -{C^s ve)'e. (2.19) 

Equation {vq). The momentum flux is 2'kRCvqYIs ^e-, the lateral force is 2'kRC Jq pdy, and 
the component of the centrifugal force which is tangential to the sphere is {dS S5 v'^/RC)sm9. 
The meridional friction force Tg dS is given by the relations 

T = ab PS v^, V = ^Jvl + v^, T0 = rsin/3, sin/3 = vq/v, Tg = aj, ps vgv, 

where Tg is the meridional component of the total shear stress, the vector r is directed opposite 
to the total flow velocity vector vgeg + v^e^, ps is the density at the bottom of the spread layer. 
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and P is the angle between the total velocity vector and the latitude. These relations follow from 
(2.16). Substituting (2.11) into the momentum equation, we obtain 

RC{T.s ve)'t = -iC^s vl)e - \c _ RCre - sinO^s vl- (2.20) 

The bottom density ps is related to the sought-for functions ve,v^p and Tg (below, we eliminate 
S5 by using the mass integral) by 

h 2Ts mpgoA 

where qwr is given by (2.8). We see that at the same surface density of the matter, the layer 
thickness, h, is inversely proportional to A. The smaller A, the higher the accuracy with which 
the radiation-pressure force offsets the difference between the gravitational and centrifugal forces. 

The subsystem (2.19), (2.20) resembles the ID hydrodynamic equations p'^. = (pu)'^, 
pu[ = —puu'^ — {pT/mp)'^. Note the "pressure" enhancement in (2.20) compared to the formula for 
an "ideal gas" Tg S5 because of the factor Geff/A (see similar remarks in Subsec. 2.3). 

Equation (v^). The rate of change in the angular momentum of the boundary-layer annulus 
is dS (S5 f(^)J RC. Writing the angular momentum flux 2itRCvq S5 VipRC and the friction torque 
about the polar axis as 

dS RC, T^ = PS \f^^r^^, 

we obtain 

R (Ss v^)[ = -{C^ S5 ve v^)'e - RC^ ps ^Jvj + v^. (2.22) 

In (2.22), in order to determine the latitudinal friction r^, we projected the total stress r onto the 
latitude, much as was done above when the meridional friction was calculated. 



2.7. Energy Transfer and Emission by the Layer 

Let us write out the energy equation (T). In this equation the formulas for the rate of change 
in the total energy of the annulus, for the enthalpy flux (the total energy flux and the power of 
lateral forces), and for the radiative losses are used. They are 

2.RCve (K + H), H = E,^^ + pdy = ^ Z^£ZI+^ 

Jo 5mp A 

qdS, K = {l/2)i:s{vl+vl). 
Adding up these expressions, we obtain the sought-for equation for the total energy 

RC {Esum + K)', = - [Cv0 {H + K% - RCq. (2.23) 
This equation allows for the energy (including radiative energy) advection and radiative cooling. 
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3. TRANSFORMATIONS AND ANALYSIS OF THE SYSTEM 

3.1. Steady-State Spread 

In the steady-state case, Eq. (2.19) is integrable. Prom the condition of mass conservation, 
we have ^ 

- M = 2TiRCveT.s (3.1) 

i.e., the matter spreads out over the two hemispheres. Since C vg S5, is constant, from the system 
(2.20), (2.22), and (2.23) we derive 

Cj:svgv'e+^c(^'^j:s] = -RCrg- sinews v^, (3.2) 

Tg = ab PS vg ^/vf+V^, 
C S5 vg (C v^y = -R T^, = ab PS ^Jvl + v'^, (3.3) 

c^,J^^^l^^.^M±l£^:^.nc,. (3.4) 

\ 2 5 rup A / 

The prime denotes differentiation with respect to 6. Because of the integral (3.1), the number of 
unknowns decreases. Eqs. (2.23) and (3.4) reflect the absence of radiative flux from the layer to 
the star. In addition, we ignore the mechanical work between the layer and the star (tv^^'^' ~ 0), 
because the star is massive (p**"^' is great even in the underlying layer) and because the velocity of 
its surface layers is low {v^"'^ <C v^; see Subsecs. 2.1 and 2.5). 

Comparison with a— friction. We disregard the friction between adjacent annuli (or 
differential elements 56) in comparison with the friction against the bottom. Indeed, in the theory 
of disk accretion (Shakura and Sunyaev 1973) with "a-friction", the turbulent stress r between the 
annuli is ahcgp (dv^p/dx), where h is the layer thickness (see Subsec. 2.1), hcg gives an estimate 
of the maximum possible hydrodynamic interaction, a is a coefficient that is small compared to 
unity, and x is the arc length along the meridian. If we take into account this interaction, instead 
of Eq. (3.3) we obtain, for example, 

-CvgJ^s {C -R C + {ajR) {C^ p c, v'^)' = 0. 

This equation contains the second-order derivative with respect to v^. The ratio of the third and 
second terms is 

{a/ab) {cs/v^) {h/RO^f. 

Neglect of the friction between adjacent annuli is justified by the smallness of {h/ROi,)"^ in 
the shallow- water approximation (2.1). Note that the term with r<^ drops out in the angular 
momentum equation for an accretion disk - there is no friction against the bottom. 

If we set C 1 in Eq. (3.4), discard Vg compared to v^, differentiate v'^, use (3.3) to eliminate 
v'lj^, and ignore radiative losses q, then we obtain {vg H)' = RTipV^. This equation together with 
Eq. (3.3) describe the transformation of rotational energy into enthalpy. 
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3.2. Natural Scales 

Let us rewrite the system (3.2)-(3.4) in dimensionless form. In units of the steUar radius 
R = 1 and the Keplerian velocity at the stellar equator = 1, it takes the form 

U^U' + ^CU'I^^-^^ =-F0-F,f, C = cose (3.5) 
UWW = -F^ + F^f, (3.6) 



^^^'^5^[~A J = - Fcf - Erad- (3.7) 

In what follows, U = vq/v^,W = v^p/v^, and T = 2Ts/mp [v^Y the dimensionless sought-for 
functions. Eq. (3.5) can be derived from (3.2) by multiplying it by vq. Below, we give the functions 
and the notation that are used to write the system. The function 

A = A(C/, W,f-e) = l-U'^ -W'^ - Gr, Gr = BCUf"^, (3.8) 

^ ^ 47r^2 a(^p(-^^)V2)4 ^ 0.65 x 1018(M/Mq)7/^ 

[see (1.2), (2.8), (2.14), and (3.1)]. For completeness, we included the centrifugal acceleration 
which is associated with the meridional motion in expression (3.8) for A. In the radiating layer, it 
is small compared to the acceleration which is associated with the latitudinal motion: U'^ <C W"^. 
In the dark layer, the centrifugal forces are insig nificant. In (3.9) and below M = W^'^ g s'^ 
The terms Fg, F^,, Fcf, and Erad which are related to the friction, the centrifugal force (cf), and 
the radiation are given by 

= F,,^^UW^ E^^^ACVf. (3.10) 

A = \—-^ -, 5 = sine, C = cose, V = VlP + W^. (3.11) 

We see that the decrease in A leads to a decrease in Fq and F^ [see also (2.20)-(2.22)]. Physically, 
this is caused by the drop in density (2.21) at the bottom due to an increase in the boundary-layer 
thickness. 

The system of equations (3.5)-(3.7) for the sought-for functions U, W, and T includes three 
parameters, LsL/Ledd^ -S, and a^. 



3.3. Transformation to an Explicit Form 



The system (3.5)-(3.7) relates the derivatives of the unknown functions U, W, and T. For this 
system to be numerically integrated, it must be resolved for the derivatives. The functions U, W, 
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and T enter into the pressure and the enthalpy which are differentiated in Eqs. (3.5) and (3.7) in a 

complicated way via the functions A, Ge//, and Gr- Since the transformation of these equations to 
a form that contains only the derivatives of the sought-for functions turns out to be cumbersome, 
it makes sense to give the final results of this transformation. We have 

U' = {erUoo - UT-eoo)M (3-12) 

W' = uj, co = i-F^ + F,f)/UW, (3.13) 

T' = (-e„iioo + Uueoo)/d, (3-14) 

where 

_ , 4 Glff - 2G,ffGr - 2U^Gr . _ 4 Ge//(Ge// + 3G.) 

uu-u ^ A2 ' ^ ~ 5 A2 ^' 

SGrUWf ismeGeffi-Geff + 2Gr)^^^ 
= 5A2 ' ''' = -5~C A^ 

_ , GAGeff + m^^ ^ _2 7G^// + 36Ge//G. - 3G,^ 
e« - ^ ^2 ' ~ 5 A2 ' 

e^ = 8^UWt eo = -4^^C/f, 

d = UmCt — UrCu, (3.15) 

^^00 = —"0 - Uyjiv - Fq - Fcf, eoo = -eo - ey,uj + Fyj - Fcf - Erad- 

Here, u^, u^,, Ut-, e^, e^, and Cr are the coefficients at the derivatives U',W', and T' in Eqs. (3.5) 
(n) and (3.7) (e), respectively. Wc represented Eqs. (3.5) and (3.7) as UuU' + UtT' = uqq and 
Bull' + Ct-T' = eoo and then solved this linear system for U' and T'. 



3.4. Phase Space and its Section by the 'Levitation' Surface (A = 0) 

The four-dimensional {6,U,W,T) phase space of the system (3.12)-(3.14) is filled with the 
integral curves U{6),W{9),T{9). Each "trajectory" (along 9) of this stream of integral curves 
represents a steady flow of the accreting plasma over the stellar surface S (Figs. 2-4) which locally 
satisfies the differential balances of mass, momentum, angular momentum, and energy (2.19), 
(2.20), (2.22), and (2.23). It is necessary to study and classify all these flows. 

An analysis shows that the numerical integration is stable only in the direction of increasing 
9. The initial data 9o and Uq = U{9o), Wq = W{9q),Tq = T{9q) at the interface between the disk 
and the spread layer are therefore required. There are many initial points in the four-dimensional 
set 9o,Uo, Wo, and Tq. Here, 9o is the angular width of the disk in the zone of its contact with 
the layer; Uq and Wq are the initial spread and rotation velocities, respectively; and Tq is the 
initial temperature. The first narrowing of the exhaustive search for 6o,Uo, Wq, and Tq is ensured 
by the fact that 6o turns out to be insignificant if ^o ^ The second narrowing is associated 
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with the separation of Wq from the set of initial data. At the disk-layer interface, W is close to 

the Kcplcrian velocity. Accordingly, we set Wq = 1 — 5 and perform three series of calculations 
with 5 = 10^^, 10"^, and 10~^' In each series, the initial points lie on the Uq,Tq (spread velocity, 
temperature at the layer bottom) plane. 

A special place in the classification analysis is occupied by the A{9, U, W,T) =0 "levitation" 
surface of the accretion layer, where A is a combination of the gravitational and centrifugal forces 
and the radiation-pressure force [see (2.8) and (3.8)]. Let us consider the intersection of this 
surface with {U,T) planes, in particular, with the {Uo,To) plane. In each section, the A = curve 
is a hyperbola, 

Uf"^ = const (3.16) 

(we assume that 9 is small and that C = cos 9 1); see (2.8) and (3.8). There are the regions 
A < and A > above and under the hyperbola (3.16). Clearly, only the region of initial data 
under the hyperbola is physically meaningful, because in this case, the total acceleration (]wr is 
directed downward and because the gravitation presses the layer against the stellar surface. 



3.5. Attracting or Limiting Surface and One-Parameter Structure 

Let us consider the "triangular" area A > that is bound by the Uq and To axes and by the 
hyperbola (3.16). This area contains the boundary (A <C 1) and distant (A Geff) subregions. 
These subregions differ greatly in area. Accordingly, the trajectory of the common position starts 
from the subregion A ~ Ggff- Consider Eqs. (3.5)-(3.7). An important point is that Uq and Tq 
are small. This is the reason why \dlnT/d9\':^l in the distant subregion and \dlnT /d9\^l in 
the boundary subregion, implying that the "hard" or steep segments of the trajectories lie in the 
subregion which is farther from the A = surface (deep in the "triangular" area) and the "soft" 
or gentle segments lie in the boundary subregion. 

Let us study the general trajectory on which A ~ Geff at the starting point. In the 
immediate vicinity of the starting point (i.e., on a steep segment), we can assume for a rough 
estimate that the basic equation is the energy equation (3.7). A changes most rapidly during the 
displacement in 6, because it contains the fourth power of T. The remaining variables {U, W, T) 
change more slowly. We disregard their change and factor them outside the diff'erentiation. Thus, 
only one unknown function remains. Assume that C ~ 1 and Gr ~ G(.ff- In the right part 
of the energy equation, we discard the small power F^f of the centrifugal force compared to the 
power Ftp which is associated with the latitudinal friction stress r<^ and radiative losses. Taking 
A 1 — W"^ — BUT^ as the unknown function instead of and rewriting the equation for this 
new unknown function, we obtain 

A^ ^ aWyf+iA/B)U ^ ^ {A/B) 
A2 ^Gefffu 4f ■ 
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We substitute the following typical values: A/B ~ 10, i7 ~ IQ-^, T ~ 10"^, W = 1 - 5, 
10~^ < 5 < 10~^, Q!b ~ 10~^. It follows from the estimate that the first term (friction) dominates 
in the numerator of the fraction before A in (3.17). For A ~ G(,ff, the first term dominates in the 
right part of (3.17). We see that A rapidly decreases on the steep segment. This decrease causes 
A to drop to such a small value that the first and second terms in the right part of (3.17) become 
equal. 



T A 




A = 



©0 9,[ 



Fig. 14. — Qualitative structure of the field of integral curves. We see that the calculations in the 
direction of increasing 6 are stable because of the attraction to the limiting surface L. The surface 

L is a separatrix - it separates the colliding streams of integral curves. The stream from the A = 
surface repulses the curves which come from the "thickness" (from the distant region; sec the text) 
to this surface. Since the ratios vq/c and ^2Ts /nip/c are small, the relative gap (for example, 
AT/r) between the L and A = surfaces is small. 

Thus, any trajectory consists of steep and gentle segments (see Fig. 14). The steep and 
gentle segments lie in the distant and boundary subregions, respectively. On the steep segment, 
the first and second of the three terms in Eq. (3.17) are important. On the gentle segment, all 
three terms are comparable. On the steep segment, A rapidly decreases with increasing 6. As one 
makes a transition from the steep segment to the gentle one, the decrease in A becomes saturated. 
The derivative |(ilnA/(i^| is small on the gentle segment, implying that the radiative losses are 
roughly offset by the frictional heating. On the steep segment, the energy balance, if by the energy 
balance we mean the equality of the right part of Eq. (3.17) to zero, is not maintained - the 
radiative losses are small compared to the heating. This implies that the left part of Eq. (3.17) 
plays an important role here. Note that the term in the left part corresponds to the advective 
transfer of radiative enthalpy [cf. Subsec. 1.3 (IV)]. The extent of the steep segment with large 
(in magnitude) derivatives is small: A^ ^ A6 = Otr — Oo, where 6tr refers to the transition 
from the steep segment to the gentle one on a given trajectory (see Fig. 14). When the complete 
system is calculated numerically, the short steep segment is traversed in small steps (a fraction of 
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A9), while the extended gentle segment is traversed with large steps (a fraction of 0^). At ~ Btr, 
the trajectory makes a sharp turn. As can be seen, the trajectories are attracted or caught in 
the boundary subregion, which lies near^ the A = surface and do not leave this subregion. For 
this reason, the "levitation" surface is called the attracting or limiting surface, hence the stability 
of the calculations and the relative "crudeness" of the results with respect to Uo and To and, 
consequently, with respect to the disk parameters^. 

Thus, the steep segments along with the distant subregion turn out to be unimportant. Only 
the boundary subregion remains. Accordingly, we can immediately take the initial data near the 
curve (3.16). In this way we further narrow down (for the third time, see above) the set of initial 
data {eo,Uo,Wo,fo). 

Let us introduce the polar (angle/radius) coordinates in the {Uo,To) plane instead of the Uq 
and Tq coordinates. Let the angle be related to the ratio 




We use the deviation 5A from the hyperbola (3.16) instead of the radius. One^ parameter i (3.18) 
now runs the set of initial data, because the deviation 5A is insignificant. A certain trajectory 
(steady accretion flow) corresponds to each value of i. We thus arrive at a one-parameter family 
of solutions, which we will call the Id family. In phase space, the integral curves from this family 
cover the surface L near the A = surface (see Fig. 14). 



3.6. Within the One-Parameter Family. The Dependence on Latitude. Critical 

Regime 

Consider an arbitrary steady flow from the family of solutions for a given Af but with 
different i. It represents the 0— distribution of U, W, and T. Let us study this flow. There are two 
characteristic regions above and below 0*. The lower region is the radiating belt. In this belt, 
[H^(^)]^ is great; the loss of rotational energy maintains the radiative flux q{9) at such a level that 
A <C 1 (or q « Qeff] see Sec. 1). The rotational energy is exhausted inside the belt. As a result, 
the luminosity q{9) abruptly decreases in the transition zone at 6 ^ 9^ - the accretion-flow surface 
darkens. Concurrently, the function A{9) increases [A{9) 1 for > 6**], g^r is large {gwr ~ 5o)) 



^ Accordingly, in this part of spread flows, the effective or total pressing is small, and the layer is thick. 

* Thus, of all the disk parameters, only the rate of mass inflow M is important for the flow in the spread layer. 
The regulation (selection) within the one-parameter (Id) family (see below) of the flow that is actually realized is 
associated with the conditions in the spreading flow at the interface between the radiating and dark layers rather 
than with the disk parameters. 

^ We obtain one value of i instead of the pair of Uo and To. The initial points lie in a narrow band near the 
hyperbola. 
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and the layer under the strong press contracts and becomes very thin outside the radiating belt 
{h is a few centimeters). 

The lower (in latitude) region is of greatest interest. Since the radiation pressure in this 
region is large, the speed of sound Cg — \/vl P is high compared to the meridional spread velocity 
VQ. For this reason, the meridional distributions inside the radiating belt become hydrostatic. 

In any solution from the family of solutions, the Mach number. 



which was calculated for a radiation-dominated plasma (Weinberg 1972), must increase with Q 
(the pressure decreases, and the flow accelerates). Indeed, the quantity (2.11) decreases with 
increasing 0, because its gradient counteracts the tangential component of the centrifugal force, 
which tends to return the plasma to the equator [see Eq. (2.20)]. 

Let us compare the Q dependences for various values of % that run the family. The Mach 
number Mag (3.19) at the interface between the disk and the spread layer [Ma6i(^? = 0)] is smaller. 



if i is larger \i is the reciprocal of the Mach number without a factor ~ JGeff/^; cf. the 



definitions (3.18) and (3.19)]. We take Mae(0o) a-t the interface between the disk and the spread 
layer and assume that the disk is thin (6*0 <C 6'^). We can thus write yia.0{6 = 9q) = Mag(0 = 0). 
At a fixed i, the function Ma^i(0) monotonically increases with 9 in the interval < 6 < 6^,. 
Consequently, Ma0{0) reaches a maximum in this interval at its edge. An analysis of the effect 
of i shows that Mae{6^), together with Ma,0{O), increases with decreasing i. On the i axis, the 
family fills the semiaxis < i < oo. The edge corresponds to the critical spread regime, in 
which the sonic point 6c lies at the boundary of the radiating belt [6c = 6^,, Mae(^^) ~ 1]. M.a0(6) 
also monotonically increases with 6 for i < i^. If i < i*, then the rotational energy at 6c, at which 
Ma0{6c) — 1, has not yet been exhausted: v^p{9c) > 0, i.e., the latitude 6c is inside the radiating 
belt. At the sonic point 9c, the solution becomes two-valued ('tips over'). Accordingly, we reject 
the solutions with z < as nonphysical. 

The critical regime occupies a special place. It restricts the Id family of subsonic (for 9 < 9i,) 
flows. In this regime, the plasma is ejected with the speed of sound from the radiating belt. 
Accordingly, the region aX 9 > 9i, has no eflFect on the flow in the radiating belt. At i > i^,, the 
regions 9 < 9-^ and 9 > 9^, are coupled. For example, the action of lateral forces (and friction tq) 
between elements 69 extends from the region 9 > 9^, into the region 6 < 6^,. 




(3.19) 




4. PHYSICS OF DECELERATION AND SPREAD 
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4.1. Numerical Simulations 

The system (3. 12)- (3. 14) was integrated numerically. In doing so, we approximated the 
derivatives by finite differences. We specified i to calculate the initial point and found the 
intersection point with the coordinates (?7oi,Tbi) of the curves i = i{U,T) and A{U,T) = 0. We 
took Uq = 0.9Uoi and Tq = 0.9Tbi as the initial values. The values of and Wq were varied 
(10-3 < ^0 < 10"^ 10-3 < 1 _ VFo < 10-1). The steep segment was traversed with the step 
~ 10-^ A^. We took 9, U, W, and T at its end as the initial values when we calculated the gentle 
segment with a step ~ 10-^^*. The round-off error was 10-^^ (double accuracy). A preliminary 
estimate of 0* was obtained from formula (1.1)' for the Eddington scale. 

The family consists of subsonic (subcritical) and supersonic (supercritical) subfamilies in the 
radiating belt. In the solutions from the supercritical subfamily, the sonic point lies inside the 
radiating belt (see Subsec. 3.6). As was stated, these solutions are discarded. 

When we say below that the fiow is subsonic, we mean that it is subsonic in the radiating 
belt. The determinant d{9) (3.15) was calculated along an integral curve. The first zero of the 
function d{9), when running in 9 from is called the sonic point 9c- At this point, the integral 
curve makes a turn, and the solution becomes two-vahicd. Since the determinant d appears in the 
denominator of the derivatives (3.12) and (3.14), the derivatives U' and T' have a singularity at 
9c- Thus, the solution exists in the interval 9q < 9 < 9c- 

4.2. One-Parameter Family as a Whole and the Passage Through the Critical Point 

As was shown in Subsecs. 3.4-3.6, the integral curves with general initial data form a 
one-parameter (or Id) family. An example of a family is given in Fig. 15 (Mag profiles, 
transcritical and supercritical solutions, 100 < i < 160) and Fig. 16 (S5 profiles in the same 
family, subcritical and supercritical solutions, 100 <i < 500). At a fixed 9, Ma^t increases, while 
Tjg decreases with decreasing i. In Fig. 15, the upper and lower solutions refer to i = 100 and 
160, respectively. In Fig. 16, the upper and lower solutions refer to i = 500 and 100, respectively. 
The supercritical solutions terminate at 9c < The subcritical solutions extend into the region 
9 > 9-1,, i.e., their termination point 9c lies to the right of 9i,- In Figs. 15 and 16, the two upper 
(Fig. 15) and two lower (Fig. 16) solutions are supercritical. The arrow in Figs. 15 and 16 marks 
the critical solution with i = and 9c = 9^,. In Fig. 15, the subcritical solutions with i > ii, run 
below (above in Fig. 16) this solution. The solutions for which i lies in the vicinity of the critical 
value (|i — ^ i^). are called transcritical. 

Let us discuss the transcritical behavior. The function 9c{i) increases with increasing i (see 
Fig. 17 and Figs. 15 and 16). This function appears to have a discontinuity at i^. In our example, 
« 151.404. At least in our example at the accuracy we chose (1.5— mm step), the discontinuity 
cannot be reduced by refining i (z_ = 151.403, i+ = 151.405). The discontinuity boundaries in R9 
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Fig. 15. — One-parameter family of steady spread flows (100 < i < 160). The dependences 
log[Mag{R9)] are for M = 4 • 10^^ g s-\ I - Wq = 10"^ 9o = 10~^, and at = W^. R9 is 
the arc length along the meridian in km, R is the stellar radius. Here i? = 10 km and M = IAMq. 
The calculations were performed using Newton's approximation. Recall that the curves terminate 
at the sonic point 9c, at which the determinant is d{9c) = [see Eq. (3.15) and Subsec. 4.1], or are 
bounded by the computational interval. 
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Fig. 16. — T,s profiles for various values of i. The four upper profiles rest on the end of the 
computational interval. The remaining profiles terminate at the sonic point. The parameters of 
the variant are given in the caption to Fig. 15. 
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are: RO. = 4.1365 and Re+ = 5.8448 km; the relative width is A9 = 2(6*+ - 9^)/{e+ + e^) = 0.34. 
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Fig. 17. — Passage of the sonic point 9c through latitude 9^ as i is varied (3.18). The latitude 9c is 
given in degrees. The parameters of the variant are given in the caption to Fig. 15. 

The latitude 9^, divides the flow into two characteristic regions (Sec. 1 and Subsec. 3.6). It 
separates the rotating (separation in W), hot (separation in T), levitating (separation in A and 
h), and radiating (separation in q) region from the nonrotating, cold, pressed and dark region. 
The separation into these regions is clearly seen on the subsonic profiles in Figs. 15 and 16 (see 
the arrows 9-^). The transition between them is sharp. We can also determine 9^, from the left 
boundary of the discontinuity in Figs. 15 and 16 (see also Fig. 17). It is marked by the arrow in 
Figs. 15-17. The discontinuity in 9c{i) is produced by the passage of the sonic point through the 
belt boundary 9^,. 

In the momentum equation (3.5), the momentum flux in the subsonic flow is small compared 
to the pressure gradient. The gradient is balanced either by the component of the centrifugal force 
Fcf or by the meridional friction Fq. Calculations show that the meridional friction^ Fg and the 
force Fcf can be disregarded inside and outside the radiating and rotating belt, respectively. This 
is yet another characteristic in which there is a clear separation at 9^,. 

Dark part of the layer. The sonic point and its position. In the dark region, the 
system (3.5)-(3.7) simplifies greatly. Here, we may assume that A = Gg// = 1 and ignore the 
rotation of matter. The flow as a whole (the radiating and dark layers) is numerically analyzed 
below. The numerical analysis must be preceded by a qualitative study. Since the variability 
of cos 9 in the equations is of little importance for a qualitative analysis, we set cos 9 = 1. The 
qualitative analysis leads us to conclude that the flow terminates at the sonic point, as confirmed 



® Of course, the latitudinal friction in the equations of angular momentum (3.3), (3.6) and energy (3.4), (3.7) 
cannot be ignored in the rotating belt. 
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below by the numerical calculations which are free from the assumed simplifications. This 

conclusion is important in constructing the ensuing theory of slow spread and settling of matter 
over a neutron-star surface, because in this way the radiative-frictional mechanism of spread 
deceleration is revealed. 

As a result of the simplifications above, we obtain 

(5[/2 - Af)U' + AUf' = -5abU^/f, (4.1) 

5UU' + 14f ' = -5A U f'^. (4.2) 
Resolving the system (4.1), (4.2) for the derivatives, we obtain 

U' = im'^ {2Af^ - 7abU) /fdi , (4.3) 

f = 5U[ {AT - 5U^)Af^ + 5abU^]/fdi, (4.4) 
di = 50U^ - 56f. (4.5) 
Eliminating 9 from the system (4.3), (4.4), we arrive at the equation 

dt ^ {Af-5U^){A/af,)f^ + 5U^ 

dU 2U[2{A/ab)f^ -7U] ^ ' ' 

in the phase plane, which allows an effective qualitative analysis, because it is two-dimensional. 

Let us estimate the ratio of the frictional (ajf/) and radiative (AT^) terms in Eq. (4.6). For 
Mi7 ~ 1, Ofe ~ 10-3, T ~ 0.5 keV, [/ ~ VT they are comparable. We offset the large (~ 10^3) 
constant A/ub by changing the scale 

U ^ \^U, f^Xf, A = {ab/Af/^ = 2.7 • 10"^ aj^^ M^/^ 

[see (3.9) and (3.11)]. Note that the velocity and temperature atU = T=l are vo = 330 Ofg^^ M^/^ 
km s-^ and T = 0.57 ag''^ M^f^ keV, where we set = lO-^ag. Replace U,f by Ma0,f. In the 
dark flow, we have Ma^ = {3/5)U'^/f [see (3.19)]. Making the replacements X = ^/W/9Mae and 
Y = T^/"^^ we can smooth out the fourth power of temperature from Stefan-Boltzmann's law, 
which complicates the analysis. In these variables, Eq. (4.6) takes the form 

dY^^Y 72/3(8 - 15X^)Y + 15X^ 

dX X 72/3(8 + 15X2)y-56X- 15X3' ^ ' ' 

A phase portrait of Eq. (4.7) is shown in Fig. 18. The function [X]c(i) = X[0c{i)] = X{i) 
gives the Mach number at the end of the integral curve (at the sonic point). It is constant 
[X{i) = 756/75] at z > z^, i.e., when the sonic point lies in the dark region. In the radiating 
region, in which z < i+ and 6c < this function slightly changes with i. For 100 < z < i^, X{i) 
is smaller than 7^6/75 by 5% and abruptly increases to this value at i*. The vertical line 4 
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0.2 0.4 0.6 0.8 1.0 



Fig. 18. — Phase plane of Eq. (4.7): 1 - the isocline of infinities, the arrows indicate the orientation 
of integral curves on the isoclines; 2 - the isocline of zeros; 3 - its asymptotic limit for y ^ oo, 4 
- the sonic point (di = 0, Mag = ^56/75). The integral curves 5 were taken from the calculations 
of the complete system which are shown in Figs. 15-17 (the subsonic subfamily, i > i*). 



{X = \/bQ/Tb) in Fig. 18 is a locus of sonic points Oc- On this line, the determinants d (3.15) and 
d\ = XT{75X^ — 56) (4.5) vanish. It is easy to show that the vertical line 4 and the isoclines of 
zeros (curve 2) and infinities (curve 1) intersect at a single point. 

It follows from the qualitative analysis of Eq. (4.7) that the integral curves which start from 
the subsonic region must necessarily intersect the vertical line 4. Indeed, the segment of the 
integral curve at 6 > 6i, starts from the subsonic region. The derivative dY/dX (4.7) is positive 
above the isocline of infinities at X < ^8/15, while the derivative dY/dO is negative [see Eq. (4.4)]. 
Hence, on the integral curve in the X, Y plane, the motion occurs in the direction of decreasing 
Y and X as 6 increases. The Y semiaxis cannot be intersected. The integral curve 5 must then 
necessarily intersect the isocline of infinities (see Fig. 18). Below the isocline of infinities, dY/dX 
is negative. This time the X semiaxis cannot be intersected. The isocline of zeros 2 lies behind the 
"triangle" that is formed by the segments of the X semiaxis, the vertical line 4, and the isocline 
of infinities 1. Therefore, the stream of integral curves 5 arrives at the vertical line 4. The curves 
crowd together near the limiting curve that corresponds to i^. The initial value of X (or Mag) 
decreases with increasing i. Thus, we established the causes of the "tipping" and determined the 
position of the sonic point in the dark region by analyzing the simplified system (4.1) and (4.2). 
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4.3. Preference of the Critical Regime. The Pedestal Under the Radiating Belt 

In Figs. 1-4 we assumed that the surface S is equipotential ("horizontal"). The equipotential 
surface of a nonrotating star is spherical. Let us find out whether the underlying surface can 
be spherical. Calculations show that the surface density T,s{9) of the plasma which constitutes 
the spreading layer is higher in the dark region (Fig. 10). In addition, the acceleration 
g^_.f j = [1 — U'^[9) — W'^{0)]gQ in a differentially rotating radiating belt is smaller than that 
in the dark region. Consequently, the weight g(,ff{0)Tis{0) of the spreading layer per unit 
area is larger in the dark region; the pressure of the spreading layer on the underlying layers 
is also higher. We ignore the motion in the underlying layer beneath the surface S. Since the 
motion in this layer is strongly subsonic, the hydrostatic contribution to the pressure is more 
important than the dynamical contribution. Because of the difference in pressure on the surface 
S from the spreading layer under the radiating region, this surface rises to form a pedestal. Let 
hus{0) = Rs{Q) — -^5(90°), where Rs is the radius of this surface. We have 

{hus)Hs = ^ ^ ^ llTus, 

pus 90 go go 

because there is no rotation under S and because the total pressing force gQ is in action. Here, 
and pb are the pressures on the surface S under the dark and radiating regions, respectively; pjjs 
is the density of the matter under S; {cs)us and Tus sue the speed of sound and the temperature 
under S; and hus and Tus are given in cm and keV, respectively. 

In addition, the radiating belt is considerably hotter with Tijs{0) = Ts{9) (see Subsec. 2.1). 
As noted above, the radiating belt heated up the underlying layers, causing the pedestal height to 
increase by {hus)T ^ {cs)s/go- 

Let us compare the height {hus)hs + {hus)T with the thickness 

Kdlnp/dy)^!" 3 A 

of the spreading layer, (see Fig. 5). We see that the thickness of the spreading flow in the radiating 
belt is much larger than the pedestal height. This is because and g^ differ. Consequently, the 
relief of the surface S is not important for the flow in the radiating belt. 

At the same time, the height hus{0) is of the order of the thickness of the spreading dark 
layer, implying that the dark flow cannot affect the flow in the radiating belt. The flow in the belt 
must therefore be transcritical. 



4.4. Spread of Dark Matter 

The dark flow is specifled by the functions vg and T5. Its segment Oi, < 6 < Oc terminates at 
the sonic point (see Subsec. 4.2). The separation Oc — increases with a decreasing initial Mach 
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number, for example, Mae(^*). In order to reach the pole {9c > 90°), it is necessary that Mag(0*) 
be sufficiently small. 

In the model (3.2)-(3.4), the surface S is impenetrable [integral (3.1)]. Above this surface, the 
newly accreted matter is transported from the source at the equator to the sink. Let us consider 
the sink. The dark layer is thin compared to the stellar radius, h -C R. It may mix with the 
underlying fluid when traversing a long path, {smix — Si,)/h ^ 1; the boundary S terminates in 
the mixing zone. 

In another case, the mixing is ineffective, and the dark flow runs up to the pole. In this case, 
we have a " bath- filling" flow with a filling "tap" near the boundary 6^. The mixing zone then is 
located not far from the tap: Smix — is of the order of (10 — 100)/i. In the "bath", the matter 
is essentially in hydrostatic equilibrium and, naturally, moves at velocities much lower than the 
speed of sound. This matter flows under the radiating belt (especially since the pressure under the 
radiating layer is lower than the pressure under the dark layer; sec Subscc. 4.3). The boundary 
S is preserved under this belt. The "bath" filling is very slow, and the flow is strongly subsonic. 
Accordingly, the weight contribution (hus)hs to the pedestal vanishes, and only the temperature 
contribution {hus)T remains. Note the possibility of subsonic circulation under S not only due to 
the viscous entrainment (gripping) by the high-velocity layer (see Subsecs. 1 and 2.1) but also 
due to the difference in the heights of equal pressure. In the Earth's atmosphere, this leads to 
trade-wind circulation which is attributable to solar heating of the equatorial zone (Palmen and 
Newton 1969; Bubnov and Golitsyn 1995). The filling causes the fluid to rise at the velocity Va or 
to settle at the same velocity. This is the radial component of the total velocity. It is easy to see 
that Va — {l/2){h/ R){ps / pus) ve is small compared to the spread velocity in the radiating belt 
- the matter traverses the belt rapidly but settles slowly. At the parameters of the problem under 
consideration, ~ 10 m s~^. Such velocity ensures a uniform increase in the surface density over 
the entire stellar surface. A particular regime of the motion of the dark fluid affects the dynamics 
of a loose belt only slightly. Accordingly, the approach (3.1)-(3.4) is acceptable for describing the 
belt all the way to the pole. Of course, flows which run up to the pole require a special analysis. 
It is necessary to introduce the parameter Va, change the boundary condition downstream, and 
add the equation for E5 to the system (3.2)- (3.4) instead of the integral (3.1). Note also that the 
allowance for the moderate penetrability of the surface S under the radiating and rotating belt 
results in certain quantitative changes, although qualitative conclusions remain valid. This can be 
seen from the results (see below) on the effect of variations in the accretion rate M on the belt 
structure. 



4.5. Effect of Friction 

Since the intensity of radiation cooling cannot exceed the Eddington limit, the area of the 
emitting surface must be large enough to remove the released heat [see (1.1), (1-1)', (1-2), (1-2)']. 
The system under study has excess energy, if the emitting area is larger than the area of the disk 
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base. The disk accretion in the disk-layer transition zone then gives way to the spread of matter 
over the stellar surface. The flow boundary ( in the plane that passes through the polar axis) turns 
90^^ in the transition zone (see Figs. 1-4 and Sec. 1). In this case, the radial flow (disk, \vr\ ^ vq) 
transforms to the meridional flow (spread, \vr \ <C vq). Because of the delay in the deceleration, a 
reservoir filled with fiuid which retains its rotation carried from the disk is formed. Its surface 
emits radiation (radiating belt). 

To assess the role of the friction model, we varied the coefficient Of, in (2.16). The cases with 
«{, = 10~^ and 10~^/3 are compared in Fig. 19. The transcritical profiles, which refer to the 
discontinuity boundaries i- and i+, are shown in this figure (see Subsec. 4.2). Let us compare the 

transcritical profiles for these values of at,- The critical value of increases roughly proportional 
to l/ofe (by a factor of 2.7). The functions Mag(9) and U{9) decrease, while Tis{0) increases oc af, 
and oc l/at, respectively. The value of 6^, increases by 6%. The functions W{9) are virtually equal, 
differing by only fractions of a percent. Consequently, the q{9) profile [see (1.2)] is invariant to ab- 
The bottom temperature Ts rises oc (l/a;,)^/^ due to the increase in S5 [see (2.14)]. 




Re, km 



Fig. 19. — Response of the function q{9) to variations in i and a^. Curves 1, 2, and 3 refer to 
moderately supersonic (1), slightly subsonic (2), and moderately subsonic (3) spread regimes for 
ah = 10~^' Curve 4 gives a subsonic transcritical profile for a^, = 10~^/3. We see that variations in 
Oft affect the q profile only slightly. 
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The degree of levitation is A c± 10~^ — 10~^, since it responds to a variation in ai, in the 
same way T5 does. The height heff (see Subsec. 4.3) is invariant to aj, variation with the same 
accuracy as W. 

As can be seen, there are parameters which change significantly as afy is varied (these include 
U, Mag, S5) and parameters which change weakly (T5, A, 6^,) or only slightly {hcff, W, q). Note 
that in view of the behavior noted above, the tangential stress (2.16) belongs to the parameters 
which respond only slightly, because ps = {4:/3)T,s /heff oc l/a^. This also follows from Eqs. 
(2.22) and (3.3) and from the invariance of W{9). The width of the discontinuity in 9c at u (A^, 
see Subsec. 4.2) depends on M and a^. 

4.6. Effect of Initial Parameters 

It is clear from the arguments which lead to formula (1.2) that the equatorial part of the 
rotating belt radiates weakly, because the normal component of the centrifugal force is large here. 
Thus, there are two bright belts at latitudes 9-^ and — The greatest brightness is reached near 
the edge of the radiating belt. The maximum is appreciably lower than unity, because where 
W :^ 1, the subtraction of the centrifugal component is significant and where W is small, there is 
nothing to radiate. It follows from the smallness of A compared to Gg// that q ^ Qeff [see Sec. 1 
for a discussion of formula (1.2)]. 

Let us analyze the effect of the initial parameters Wq and 9o, associated with the disk, on the 
results. We compare the cases with 1 — Wq = 10~^ and 10~^- They correspond to a change in the 
deviation of the rotation velocity from the Keplerian velocity by an order of magnitude. For this 
change in Wq, the kinetic energy increases by 2%, increases by 42%, the belt width 9^, increases 
by 3%, MaQ inside the belt decreases by 12%, the surface density T,s increases by 7%, the rotation 
velocity W increases by 1.4%, the temperature rises by 0.8%, the degree of levitation A decreases 
by 13%, the height heff increases by 14%, and q decreases by 3%. 

Let us increase the meridional extent of the disk base by a factor of 10: ^0 = 10~^ is replaced 
by 10~^- The critical value of increases by 8%. The distributions of Mag,S5,T^ and other 
variables for disks of thickness 2R9o = 20 and 200 m differ only slightly. The distribution in 
9 seems to shift by an amount of the order of the difference between the values of ^o- Thus, 
for 1 — Wq <^ 1 and 9q <^ 9^,, variations in Wq and have a marginal effect on the results, in 
particular, on the q{9) profile. 

4.7. Dependence on Accretion Rate 

Above, we discussed the main points regarding variations in a;,, VFq, and 9q. Let us now 
study the effect of M. The critical value i*(M) and the belt width 9i,{M) are plotted versus the 



-45- 



accretion rate in Fig. 20. In the calculations, we approached the critical value by compressing 
the transcritical vicinity by the fork method (see Subsecs. 3.6, 4.2, 4.3). The width 6'^(M) was 
determined from the last calculated supercritical profile. We took the position of the sonic point 
in this profile as the belt width. An analysis of the effect of variations in M on the spread-flow 
structure shows that the relative discontinuity width (Subsec. 4.2) decreases with increasing 
M. At M Mpoie, the discontinuity disappears, and the dependence on i becomes smooth; the 
edge of the radiating belt comes close to the pole. Under these conditions, the edge of the rotating 
belt was taken as 9^{M). We chose the smallest value of i at which this edge lies at a distance 
from the pole and took it as i*. 




Fig. 20. — Critical value of i*, radiating-belt width and surface density {J^s)-k at the belt 
boundary versus M. 

In Fig. 20, we also plotted ij^s)* versus M. The value of {J^s)-k is Tis{Q-k)i where Oi, was 
determined from the last calculated supercritical profile. Near the point of stoppage at the pole 

for M ~ Mpgie, we took the edge of the rotating belt as 9^. The function {Tis)-k{M) begins to 
rapidly increase near the pole. It is physically clear that after the point of stoppage of the loose 
radiating belt at the pole, the spread of matter slows down greatly (Subsec. 4.4). The slowdown 
causes the surface density to increase proportionally. The power law fitting the dependences 
are ^ (S5)* [g cm-=^] ^ 35MiV^, Re^ = [km] ^ lAM^j^. 

Let us define the contrast as the ratio of E5 in the central zone to (Sg)^. The contrast 
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increases with increasing M. The fact that it is greater than unity is attributable to the existence of 

a strongly subsonic central zone. For this reason, the spreading layer stretches when it approaches 
the edge of the radiating belt. While at large M there are a subsonic center and a transonic edge, 
at small M the entire radiating belt is transonic. 

4.8. Limiting Magnetic Field 

Let us estimate the magnetic-field strength that still has no effect on the pattern of 
deceleration and spread outlined above. Compare the magnetic pressure H'^/Stt and the pressure 
PS = 9eff ^5 (2-2) at the bottom of the spread layer. To obtain an estimate, we assume that 
9eff ~ 9o- We take C^s)* as a typical value of S5 and use the fitting dependence on M which we 
derived in Subscc. 4.7. As a result, wc obtain H < Hmax, H-max = 4 x 10^(Mi7)'^"''^ gauss. A 
weaker field appears to "get buried" under the layers of accreting plasma. X-ray bursters give an 
example of neutron stars with a weak magnetic field. 

CONCLUSION 

The spread of matter during disk accretion was considered. We found a radiating belt to 
be formed in a certain range of M. The belt width Oi, depends on M. At LsL/Lgdd — 3 x 10"^, 
the width is of the order of the thin-disk thickness Oq = 10~^ 10~^. The belt disappears at 
Lsh/Ledd ~ Lpoie/Ledd ~ 0.9, becausc the entire stellar surface radiates. The belt rotates. Its 
rotation velocity depends on latitude. This belt levitates in the sense that the pressing acceleration 
Qwr is small compared to the gravitational acceleration qq. A delay in the deceleration of rotation 
is responsible for the appearance of the belt. This delay in the deceleration is attributable to a 
great decrease in the bottom density due to the levitation and loosening of the matter. 
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